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ABSTRACT 

We present a detailed comparison of the CO (3-2) emitting molecular gas between a local sample 
of luminous infrared galaxies (U/LIRGs) and a high redshift sample that comprises submm selected 
galaxies (SMGs), quasars, and Lyman Break Galaxies (LBGs). The U/LIRG sample consists of our 
recent CO (3-2) survey using the Submillimeter Array while the CO (3-2) data for the high redshift 
population are obtained from the literature. We find that the Lco(3-2) ^^"^ LpiR relation is correlated 
over five orders of magnitude, which suggests that the molecular gas traced in CO (3-2) emission is 
a robust tracer of dusty star formation activity. The near unity slope of 0.93 ± 0.03 obtained from 
a fit to this relation suggests that the star formation efficiency is constant to within a factor of two 
across different types of galaxies residing in vastly different epochs. The CO (3-2) size measurements 
suggest that the molecular gas disks in local U/LIRGs (0.3 - 3.1 kpc) are much more compact than 
the SMGs (3 - 16 kpc), and that the size scales of SMGs are comparable to the nuclear separation (5 
- 40 kpc) of the widely separated nuclei of U/LIRGs in our sample. We argue from these results that 
the SMGs studied here are predominantly intermediate stage mergers, and that the wider line- widths 
arise from the violent merger of two massive gas-rich galaxies taking place deep in a massive halo 
potential. 

Subject headings: galaxies: formation, galaxies: starburst, cosmology: observations, galaxies: high 
redshift, submillimeter 



1. INTRODUCTION 

Ultra/Luminous Infrared Galaxies (U/LIRGs) are 
sources that emit large amounts of flux in the far- 
infrared (FIR) bands, with FIR luminosities in the range 
10"-i2 Lq (for LIRGs) or even larger (> lO^^ Lq for 
ULIRGs) (see reviews by Sanders & Mirabel 1996; Lons- 
dale, Farrah & Smith 2006). These galaxies were discov- 
ered by the Infrared Astronomical Satellite (IRAS) to 
be fairly rare in the local universe, but show significant 
increase as a function of redshift (Hacking, Houck & Con- 
don 1987). FoUowup optical studies of nearby U/LIRGs 
have revealed that the morphologies of a significant frac- 
tion of them resemble interacting/merging systems (Ar- 
mus, Hcckman & Miley 1987), mostly powered by star- 
bursts with increasing contribution from Active Galactic 
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Nuclei (AGN) for IR-brighter galaxies (Kim, Veilleux, 
& Sanders 1998). These galaxies harbor large amounts 
of molecular gas (Sanders et al. 1986) that is often con- 
centrated near the compact nuclear regions (Downes & 
Solomon 1998; Bryant & Scoville 1999). These obser- 
vational findings were analyzed in conjunction with nu- 
merical simulations that date back as far as the 1970's 
(Toomre & Toomre 1972), further advanced with the in- 
clusion of gas dynamics and star formation recipes in 
the 1990's (Barnes k Hernquist 1996; Mihos & Hcrn- 
quist 1996) and with realistic feedback mechanisms im- 
plemented in the 2000's (Springel, Yoshida & White 
2001; Springel, Di Matteo & Hernquist 2005; Cox et al. 
2006). As a result, it is widely believed that tidally trig- 
gered gas compression and subsequent starbursts (and/or 
AGN fueling) heat the surrounding dust, are reprocessed 
into far-infrared emission, and result in the intriguing 
U/LIRG phenomenon we observe in the local as well as 
the high redshift universe. 

Since the late 1990's, with the advent of sensitive 
mm/submm bolometer cameras mounted on single dish 
telescopes, observations of the blank sky have revealed 
the ubiquitous presence of unresolved sources that ap- 
pear to account for a substantial fraction of the total 
infrared background radiation (see Lagache, Puget, & 
Dole 2005, for a review). By using the correlation be- 
tween radio and FIR (Condon 1992), interferometric ra- 
dio observations toward these sources and subsequent 
foUowup optical spectroscopy have provided substantial 
evidence that the origin of this emission is dusty star- 
forming galaxies that reside at z > 2 (Chapman et al. 
2004b). The FIR luminosities and molecular gas prop- 
erties in these submillimeter selected galaxies (SMGs) 
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are generally an order of magnitude larger than the local 
U/LIRGs (Grove ct al. 2005; Tacconi et al. 2006, 2008). 
While the FIR luminosities in local U/LIRGs appear to 
be powered primarily by starbursts and in some cases 
a central AGN (e.g. Genzel et al. 1998), recent studies 
have suggested that the AGN contribution to the FIR 
luminosity is negligible in the high redshift SMGs (Pope 
et al. 2006, 2008). 

It has long been proposed that U/LIRGs and quasars 
are linked through dynamical evolution of the host sys- 
tems, likely by a merger event (Sanders et al. 1988a; 
Hopkins ct al. 2005). The initial stages of the merger 
could be dominated by dusty starbursts (i.e. U/LIRGs 
and SMGs), until the central nuclear activity is stimu- 
lated by massive gas accretion. When the cold gas is 
consumed either by starbursts or accretion onto the cen- 
tral AGN, the quasar may become visible as interven- 
ing dusty high column density material along our line of 
sight clears away. This merger evolution scenario is par- 
tially supported by recent deep and high-resolution imag- 
ing of quasar host galaxies (Dunlop et al. 2003; Sanchez 
et al. 2004; Veilleux et al. 2006; Zakamska et al. 2006), 
where morphological signatures of interactions/mergers 
and spheroidal systems are evident. 

Since the important fuel of the activity, whether it be 
starbursts or AGN, is warm and dense gas, a large sample 
of molecular line data in IR bright galaxies and quasars 
at different epochs is needed. This data set allows us 
to test whether the evolutionary scenario is consistent 
with observations in the local universe, and then to test 
whether it holds true at high redshifts. Single dish and 
interferometric observations of the CO (1-0) emission in 
local U/LIRGs have been carried out previously, provid- 
ing significant evidence that the amount of gas in these 
systems is massive, and highly concentrated in the nu- 
cleus (e.g. Downes & Solomon 1998). Similar observa- 
tions have been conducted in CO emission of the red- 
shifted high- J transition lines from SMGs/quasars (see 
Solomon & vanden Bout 2005, for a review), but these 
observations are difficult due to the requirement of a pre- 
cise knowledge of the redshift and the lack of sensitivity 
at mm/submm bands, and therefore the sample is biased 
toward the brightest (or gravitationally lensed) sources. 

The majority of the molecular gas detections in the 
bright SMGs/quasars thus far from z > 2 sources are 
CO (3-2) (or higher J) line redshifted to the mm bands. 
However, directly comparing the properties of the CO (3- 
2) emission in the SMGs/quasars with the CO (1-0) 
emission in the local galaxies can be significantly bi- 
ased by excitation effects^^. Thus, to properly assess 
the characteristics of gas in U/LIRGs and high redshift 
SMGs/quasars, a thorough comparison of a single J tran- 
sition emission line is needed. We have carried out a pro- 
gram at the Submillimeter Array (SMA)^^(Ho, Moran & 
Lo 2004) to observe 14 U/LIRGs in the local universe in 
the CO (3-2) and CO (2-1) lines (Wilson et al. 2008, Pa- 
per I hereafter). Because the critical density of CO (3-2) 
emission is relatively high, it is considered as an impor- 

11 Tex = 33 K and nc„t = 3.6 X 10"' cm-^ for CO (3-2) whereas 
Tex = GK and iierzt = 2.2 X 10'^ cm-'^ for CO (1-0) 

i'^ The SubmiUimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



tant tracer of the physical condition of the star forming 
gas and the associated kinematics. This common belief, 
however, has not been fully tested observationally be- 
cause of the lack of high resolution observations of nearby 
and distant sources. Recent high resolution studies of lo- 
cal U/LIRGs show, in some cases, that the distribution 
and the peak of the CO (3 2) emission could be different 
from the lower transition CO (e.g. lono et al. 2004). 

In order to conduct an analysis that is not biased by 
excitation effects, we present a detailed comparison of the 
CO (3-2) emitting molecular gas between a local sample 
of U/LIRGs and a high redshift sample that comprises 
SMGs, quasars, and two Lyman Break Galaxies (LBGs) 
that have been recently detected in CO (3-2). We also 
include CO (3-2) data from a sample of local quiescent 
galaxies in some of the analyses. The U/LIRG sample 
consists of our recent CO (3 2) survey using the SMA 
(Paper I) while the CO (3-2) data for the high redshift 
and local quiescent galaxies are obtained from the liter- 
ature. 

In §2, we present the sample data and the formalism for 
deriving the relevant parameters. In §3, we investigate 
the correlation between CO (3-2) and FIR luminosities in 
a variety of sources, and offer possible interpretations of 
the derived slope and the scatter seen in the correlation. 
We compare the CO (3-2) source size in different sources 
and offer clues to the merger stage of SMGs and quasars 
in §4. We present our conclusions in §5. We adopt Hq 
= 70 km s-i Mpc-\ Qm = 0.3, Oa = 0.7 for all of the 
analysis throughout this paper. 

2. SAMPLE DATA AND DERIVED QUANTITIES 

Observational details and the properties of the local 
U/LIRG data are described in Paper I. The U/LIRG 
sample consists of 14 sources, of which six are interact- 
ing/merging galaxy pairs whose nuclei are resolved into 
multiple components ( "pre-coalesccncc" mergers), and 
eight are mergers whose nuclei are unresolved even with 
space- based optical images ("post-coalescence" mergers). 
For the remainder of the paper, the local infrared bright 
galaxies with Lfir > 10^^ L© are referred to as ULIRGs, 
whereas galaxies with LpiR < 10"'^^ L0 are referred to as 
LIRGs. The term U/LIRG is used when referring to the 
entire sample presented in Paper I, regardless of the FIR 
luminosity. 

Interferometric CO (3 2) data toward SMGs were ob- 
tained from published results in Neri et al. (2003), Sheth 
et al. (2004), Greve et al. (2005), Tacconi et al. (2006, 
2008) and Knudsen et al. (2008). A few sources that 
were observed in Greve et al. (2005) were also observed 
by Tacconi et al. (2006) at higher angular resolution, 
and we adopt the latter for our analysis. In addition, 
CO (3 2) emission in MIPS-J1428 (lono et al. 2006a), 
a submm bright HyLIRG (Hyper Luminous InfraRed 
Galaxy; Lfir > 10^^ L©) at 2: = 1.3, was also included 
in our SMG sample. The CO (3-2) data for the quasars 
were compiled from Hainline et al. (2004), Walter et al. 
(2004), Solomon & vanden Bout (2005), and Coppin et 
al. (2008). The two CO (3 2) detected LBGs, MS 1512- 
cB58 (cB58 hereafter) and J213512-010143 (a.k.a. the 
"cosmic eye", J213512 hereafter), were obtained from 
Baker et al. (2004) and Coppin et al. (2007), respectively 

The far infrared luminosities (Lfir) are obtained from 
the IRAS Revised Bright Galaxy Sample (BGS) (Sanders 
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ct al. 2003) for all of the U/LIRGs. The LpiR of 
SMGs/quasars arc derived by fitting a theoretical star- 
burst spectral energy distribution (SED) (Efstathiou, 
Rowan-Robinson, & Siebenmorgen 2000) to the observed 
FIR to radio flux densities (Yun & Carilli 2002). The 
same SEDs were used to compute the Lpm of the three 
local ULIRGs in our sample, and wc found that they 
are consistent with the Sanders ct al. (2003) values to 
within < 50%. The LpiR of the LBGs are obtained from 
Baker et al. (2004) and Coppin et al. (2007). The sam- 
ple sources along with various physical properties used 
in this study arc presented in Table 1. 

We derive the following physical properties for the 
analysis of this paper. The derived physical properties 
of all of the sources are presented in Table 2. 

(1) The CO (3-2) luminosity, ^'^0(3-2) [K km s^^ pc^], 
is derived using 

4o(3-2) = 3.25 X 10^ ScO(3-2) i^^ts (1 + ^)"' Dl (1) 
where 500(3-2) is the integrated CO(3-2) intensity in 
Jy km s~^, fobs is the observed frequency in GHz, Dl is 
the luminosity distance in Mpc, and z is the redshift of 
the source (Solomon & vanden Bout 2005). 

(2) The full width at half maximum (FWHM) of the 
CO (3-2) line are estimated using the integrated spectra 
shown in Fig. 25-29 of Paper I, and obtained from the 
published results for the SMGs, quasars, and LBGs. 

(3) The source sizes of U/LIRGs arc the deconvolved 
FWHM diameters derived by fitting a 2 dimensional 
Gaussian to the CO (3-2) integrated intensity maps. 
When the deconvolution failed in some of the sources 
(mainly due to low S/N), the major and minor FWHM 
axes obtained from the Gaussian fits in the image do- 
main are given as upper limits to the source size. All 
of the U/LIRGs are resolved with the SMA beam, ex- 
cept for IRAS 17208-0014 whose deconvolution produces 
a point source. However, visual inspection of its inte- 
grated intensity map in Figure 2 of Paper I suggests that 
the CO (3-2) emission appears to have a bright compact 
component, as well as a resolved extended structure. 

The published sizes are used for the 6 SMGs that 
are resolved by the interferometer. Most of the 
quasars are unresolved except for SDSS J1148-I-5251, 
IRAS F10214-I-4724, and Cloverleaf. High angular reso- 
lution VLA observations (Walter et al. 2004) have suc- 
cessfully resolved the structure in the z = 6.4 quasar 
SDSS J1148-I-5251, whereas OVRO observations (Yun et 
al. 1997) in conjunction with a lensing model have in- 
ferred a source size of 0'.'3 for the highly lensed z = 2.6 
quasar Cloverleaf. For the remaining unresolved SMGs, 
quasars, and LBGs, the minor axis of the beam (FWHM) 
is used as an upper limit to the source size. 

Because the physical resolution of the high redshift 
population is worse than the local LIRGs, some of the 
widely separated sources that are resolved with the SMA 
may not be resolved using 0'.'8 resolution when they are 
placed at ^ = 2. We have tested this by convolving the 
CO (3 2) images of all of the widely separated pairs with 
O'.'S (^ 8 kpc at z = 2) resolution, and found that all the 
galaxies are clearly resolved into multiple components 
except for Arp 299, which is only marginally resolved. 

(4) The star formation eflaciency (SFE), 

Lfir/L(-,qj3_2), is derived by taking the ratio be- 



tween FIR luminosity and CO (3 2) luminosity. The 
ratio between FIR luminosity and molecular gas mass 
(i.e. Lfir/Mhj) is often used to infer the SFE of 
a galaxy, but we opted to use the luminosity ratio 
(Lfir/Lqq^3_2)) this study to avoid introducing ad- 
ditional ambiguities through variation and uncertainties 
pertaining to the conversion from Lf.Q^3_2) to Mhj- 

We caution that the CO (3-2) properties at low and 
high redshifts are generally lower limits for the following 
reasons. In synthesis imaging, incomplete uv coverage 
results in decreasing sensitivity to extended structure. 
While the SMA observations are sensitive to structure 
with size scales of < 12" — 16" (7 10 kpc on average for 
these galaxies), the missing flux calculated by comparing 
to single dish data is 2 - 66% (Paper I). Thus the diffuse 
and extended CO (3-2) emission would be undetected in 
the SMA observations, and the values derived in this pa- 
per are lower limits to the physical quantities intrinsic to 
the source. On the other hand, observations of high red- 
shift galaxies are limited in both sensitivity and angular 
resolution. The relatively coarse angular (physical) reso- 
lution of < ^.'8 (^ 6 kpc) allows us to analyze images that 
do not suffer significantly from the missing-flux problem. 
However, the high redshift population suffers from lim- 
ited sensitivity (with mass sensitivity of ~ 10^" M©), 
and any extended CO (3-2) emission could be buried in 
the noise. Therefore, the physical quantities derived here 
from CO (3 2) data for both the low and high redshift 
populations are likely lower limits. 

3. THE CO (3-2) LUMINOSITY AND THE STAR 
FORMATION EFFICIENCY 

The luminosity of the CO (3-2) emission line can be a 
direct measure of the amount of dense molecular gas fuel- 
ing nuclear (or extended) star formation and/or the cen- 
tral AGN. The average CO (3-2) luminosities (Lqq(^-^_2)) 
derived in U/LIRGs, SMGs, and quasars are (2.6±0.5) x 
10^, (4.4±l.l)xl0i°and (5.0±1.0) x 10^° K km s'^ pc^, 
respectively. The L^Q^^ g) fo'" the two LBGs are 4.4 x 10^ 
and 2.9 x 10^ K km s'^ pc^ for cB58 and J213512, respec- 
tively. The average CO (3-2) luminosity in the U/LIRGs 
is more than an order of magnitude lower than in the 
SMGs and quasars, but comparable to the two LBGs. 
Under the assumption that the powering source of the 
large Lfir is mostly from starbursts (see §3.1), the com- 
parison between ^qq(^^_2) Lfir relates the amount 
of available dense molecular gas to the amount of current 
massive star formation traced in FIR dust emission. This 
analysis assumes an environment where the gas is suffl- 
ciently dense and warm to thermally (or radiatively) pop- 
ulate the J=3 rotational energy level of carbon monoxide, 
and that the gas and dust are coupled spatially. We in- 
vestigate the relation between ^co(^3_2) Lfir in the 
following subsection. 

3.1. The -^co(3-2) ^ ^FIR, Relation 

In Figure 1 {left), we compare the CO (3-2) line lumi- 
nosity with the FIR luminosity for all four sample pop- 
ulations. In addition, we have plotted the same lumi- 
nosity relation for a sample of local galaxies published in 
Mauersberger et al. (1999) and Komugi et al. (2007) to 
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increase the luminosity range. The L(~,q^3_2) " Lfir re- 
lation including all four populations (and local galaxies) 
is correlated over five orders of magnitude in luminosity, 
which suggests that the molecular gas seen in CO (3-2) 
emission is a robust tracer of star formation activity. A 
least-squares fit of the form log Lco(3-2) ~ *^ Lfir+/3 
gives a = (0.93 ±0.03) and (3 = (-1.50 ± 0.33). Exclud- 
ing the local galaxies or the quasars, which may have 
significant dust heating due to the central AGN, does 
not change the slope with significance. One potential 
concern here is that the large range in distance may 
artificially produce a correlation in log-log space. To 
check whether the correlation seen here is statistically ro- 
bust, we employed the partial Kendall r-cocfficicnt with 
censored data (Akritas & Siebert 1996) where wc used 
the luminosity distance as a test variable. This analysis 
yields a probability of 7.3 X 10~^ that these two vari- 
ables will produce a false correlation after the effect of 
distance is removed, providing solid statistical evidence 
of correlation in logarithmic space. 

The near unity slope derived from the Lfir, - ^co('i-2) 
relation suggests that the efficiency of converting CO (3- 
2) emitting molecular gas to massive stars (i.e. SFE) is, 
within a factor of two (i.e. from the standard deviation; 
see below) , nearly uniform across different types of galax- 
ies residing in vastly different epochs. To illustrate this 
better, we plot the SFEs against the FIR luminosities in 
Figure 1 {right). The average SFEs in the various galaxy 
samples range from a low of 170 for the SMGs and LBGs 
to 250 for the LIRGs to 430 for the quasars and a high of 
580 for the ULIRGs. The combined average of the SMGs, 
quasars, LBGs and U/LIRGs is 375 ± 42 (standard de- 
viation = 327) L0 (K km s^^ pc^)^^. The average 
SFE of the SMGs is lower than the U/LIRGs by a factor 
of two, while a slightly higher average SFE seen in the 
quasars (and also the ULIRGs) is possibly attributable 
to higher dust temperature as a result of dust heating by 
the central AGN. Interestingly, the SFEs of the SMGs 
are more consistent with the widely separated LIRGs 
(and the LBGs and local galaxies) than the ULIRGs and 
quasars. We will return to this point in our discussion of 
the evolutionary status of the SMGs in §4.2. The SFE 
and LpiR, properties of the LBGs are also generally simi- 
lar to the LIRGs, but this sample is currently limited to 
only two galaxies. 

In U/LIRGs, the exact fraction of the AGN/starburst 
contribution to the FIR luminosity appears to depend 
on the source (Armus et al. 2007). Half of the four- 
teen U/LIRGs in our local sample show evidence for an 
AGN, but only in the two ULIRGs Mrk231 and Mrk273 
is it possible that the AGN makes a significant contribu- 
tion to the bolometric luminosity (Paper I, and references 
therein). This variation suggests that the scatter in the 
SFEs of U/LIRGs not only reflects the physical char- 
acteristics of each galaxy, but also the AGN heating of 
Lfir, in some sources. However, we mention the overall 
trend that the FIR output in local IRAS sources is dom- 
inated by starbursts and not AGN (Yun et al. 2001), as 
is evident from the robustness of the empirical radio- FIR 
correlation. In addition, while there is evidence from X- 
ray surveys that a good fraction of SMGs contain AGNs 
(~ 80%; Alexander et al. 2005), the energetic contribu- 
tion of the AGN to the FIR output in most of the SMGs 



appears to be minor (Pope et al. 2008; Valiante et al. 
2007). Indeed, there is even evidence that most of the 
far-infrared luminosity in quasars comes from a central 
starburst (Lutz et al. 2007; Wang et al. 2007; Coppin 
et al. 2008), although there are exceptions (Weiss et al. 
2007). 

In Paper I, it was found that the relation between 
Mh2 (or L(^Q^3_2)) Lfir was not correlated with 
any significance, leading to an argument that the degree 
of gas concentration determines the level of star forma- 
tion activity rather than the overall amount of molecu- 
lar gas mass available in the galaxy. Due to large scat- 
ter in the SFEs of LIRGs (see Figure 1 {right)), the 
correlation indeed breaks down if we only consider the 
U/LIRGs in the fit. The LIRGs that have particularly 
low SFEs (< 150) are NGC 6240, Arp 193, VV 114, 
NGC 5257/8, and NGC 5331. Optically thin CO (3 2) 
is suggested from previous analysis on NGC 6240 (lono 
et al. 2007), and therefore the CO (3-2) emission may 
be over-luminous with respect to other LIRGs with op- 
tically thick CO (3-2) emission. Although Arp 193 is 
apparently a late stage merger, W 114, NGC 5257/8, 
and NGC 5331 show large nuclear separation suggestive 
of early stages of interaction. The SFEs of these galaxies 
are 93 - 117 L0 (K km s~^ pc^)^-'^, which is a factor of 
three lower than the average of the entire sample. The 
factor of two lower average SFE seen in LIRGs (which 
have a mix of compact and widely separated sources) 
compared to ULIRGs (which are mainly compact) may 
suggest that the global SFE is lower in early to interme- 
diate stage pre-coalescence mergers than in the centers 
of more advanced stage mergers, likely because the gas 
in intermediate stage mergers has not settled in a steady 
state in the rapidly evolving host stellar potential. 

Finally, we note that the average SFE in lo- 
cal U/LIRGs is a factor of seven larger than 
Lfir/L(-,q^2-i) ^ 50 found in two z ~ 1.5 BzK- 
color selected ULIRGs (Daddi et al. 2008). How- 
ever, Lfir/L(-,q|.2_j-) may be significantly lower than 
Lfir/Lqq^3_2') in the same galaxy owing to excitation 

effects, and future CO (3-2) observations of the same 
z ^ 1.5 BzK selected ULIRGs are necessary to compare 
the true CO (3 2) properties of molecular gas in ULIRGs 
at this redshift. Tacconi et al. (2008) have obtained up- 
per limits to the CO (3-2) emission in three optically 
selected star forming galaxies similar to the BzK sam- 
ple; however, the average lower limit to their SFE (> 90, 
where wc have estimated Lfir from their measured star 
formation rates using the relation in Kennicutt 1998) im- 
plies that the local U/LIRGs have at most a factor of 
three times larger SFE. Upper limits to the CO (3-2) 
emission in 5 SMGs and 3 quasars (Greve et al. 2005; 
Coppin et al. 2008) give lower limits to their SFEs of 
> 160 to > 1000. While most of these limits lie above 
the average SFE shown in Figure 1 {right), significantly 
more sensitive CO (3-2) observations of these sources 
would be needed to show conclusively that they deviate 
significantly in their properties from the detected SMGs 
and quasars. 

3.2. Slope in the Lfir - Lqq Relation 



5 



The CO (3-2) derived slope (aco(3-2) = 0.93 ± 0.03) 

in the Lfir - L^q relation is significantly steeper than 
a previous compilation which includes different excita- 
tion lines from local to high redshift sources (i.e. aco = 
0.62 ± 0.08; Grcvc ct al. 2005). This is expected be- 
cause the molecular ISM in most of the SMGs is sub- 
thermally excited in CO transitions higher than J — 4-3 
(Weiss et al. 2007b), and their sample of local U/LIRGs 
mainly consists of observations of J = 1-0. A study 
that investigates only the CO (3 2) lino in local galax- 
ies with FIR luminosities in the range 10^ -10^^ L© 
was conducted by Yao ct al. (2003), where they find 
aco(3-2) = 0.70±0.04. A more recent study (Narayanan 
et al. 2005) that relates the CO (3-2) to IR luminosity in 
a similar luminosity range suggests a near linear correla- 
tion (aco(3-2) = 1-08 ± 0.07). While there is variation 
in the slope among different studies, the slope between 
the CO (1-0) line and the FIR luminosity is even flat- 
ter (aco(i-o) = 0.58 ± 0.07; Yao et al. 2003) than that 
derived using the CO (3-2) emission alone. The rela- 
tively flat CO (1-0) slope suggests that the CO (1-0) de- 
rived SFE increases with increasing FIR luminosity (i.e.; 
Sanders et al. 1988a), and the molecular gas traced in 
CO (1-0) emission is not a linear tracer of star forma- 
tion activity. 

It has been demonstrated that the HCN emission, 
which is a higher density gas tracer than the commonly 
used CO (1-0) emission, produces a linear correlation 
with FIR luminosity (q!hcn(i-o) = 1-00 ± 0.05; Gao & 
Solomon 2004). Similar recent studies (Gao et al. 2007; 
Riechers et al. 2007; Gracia-Carpio et al. 2008) that in- 
clude HCN observations from high redshift galaxies have 
shown, however, that Lfir/Ljjqj^ in the most distant 
galaxies is larger than the local sample by at least a fac- 
tor of two, and the resultant flt between Lhcn(i-o) ^^"^ 
Lfir becomes flatter (q:hcn(i-o) = 0.81 ± 0.06; Gracia- 
Carpio et al. 2008) when high redshift galaxies are in- 
cluded. These new HCN results suggest that the dense 
gas fraction in bright high redshift galaxies is even higher 
than in the local IR bright galaxies. 

In Figure 2, wc show the different slopes (amoi) plotted 
against the critical densities of the CO (1-0), CO (3-2) 
and HCN (1-0) transitions. We calculate the critical 
densities from Auij~i^ where A^i is the Einstein A coef- 
ficient for the upper (u) to the lower (/) energy state and 
7„( is the corresponding collision rate obtained from the 
compilation of Schoier et al. (2005). A general trend that 
(Xmoi approaches unity for higher critical density tracers 
is evident from Figure 2. 

The slope derived from CO (1-0) emission, q:co(i-o)) 
could be 2/3 if a fixed fraction of gas converts to stars 
each free-fall time (Kennicutt 1998). From simple the- 
oretical considerations, Krumholz & Thompson (2007) 
argue that, because the critical density of the CO (1- 
0) line is lower than the higher density tracers such as 
HCN, the CO (1-0) line traces the molecular gas that 
has densities close to the median density of a galaxy. 
On the other hand, molecules with high critical densities 
(such as the HCN (1-0) hue) trace the ISM with much 
higher densities. In this case, LpiR/Lygj^j-^ ^,^ does not 
depend strongly on the mean density of the galaxy, and 
therefore the slope derived from the HCN (1-0) measure- 



ments (i.e. aHCN(i-o)) should be close to unity. They 
find, however, that the dependence of Lfir/Ljj(-,j^^^_qj 
to the mean density becomes more significant at higher 
(10''^"' cm^'') mean densities, suggesting that the slope 
can be flatter at the highest luminosity regime. 

We plot the two theoretical predictions for 
Q^co(i-o) (dotted line) and Q;hcn(i-o) (short dashed 
line) in Figure 2 as horizontal lines. The observations 
for these two emission lines are consistent with the 
theoretical predictions, and the new measurement by 
Gracia-Carpio et al. (2008) is also consistent with the 
prediction by Krumholz & Thompson (2007) that q:hcn 
is less than unity when galaxies with higher mean den- 
sities (rcprcscmted by high redshift sources) are included 
in the fit. Our CO (3 2) measurements, as well as those 
measured by others, fall between the CO (1-0) and 
HCN (1-0) derived slopes, but the work by Krumholz 
& Thompson (2007) was limited to an isothermal case, 
and therefore could not investigate higher transition 
lines that are excited at higher temperatures such as the 
CO (3-2) line. 

3.3. Molecular Excitation and Large Velocity Gradient 

Modeling 

A theoretical investigation was performed by 
Narayanan et al. (2007b) to understand the dynamically 
evolving ISM of isolated and interacting galaxies. Prom 
their N-body/SPH simulations that include radiative 
transfer, they find that the total CO (3-2) luminosity 
has a mix of contribution from the dense cores (sites of 
star formation) and the more extended molecular gas 
that is subthermally excited. They argue that CO (3-2) 
and FIR luminosities are linearly correlated under the 
assumption that the gas density and SFR density are 
related by a Schmidt law (Schmidt 1959) index of 1.5. 
Further, a notable scatter is seen in the low luminosity 
end of their CO (3-2) - FIR luminosity relation, which 
is caused by line trapping of the CO (3-2) line by 
molecular gas surrounding the dense star forming cores 
and suggests that subthermally excited CO (3-2) gas 
can provide a substantial contribution to the overall 

^CO(3-2)- 

In order to investigate the CO (3 2) excitation at 
different spatial scales, we derived the total {Rtotai) 
and peak [Rpeak) CO (3-2) to CO (1-0) ratios in 
U/LIRGs using the CO (1-0) data available in the lit- 
erature. The results are presented in Table 3. The av- 
erage Rtotai is 0.48, ranging from 0.20 (Arp 55) to 1.25 
(Arp 299W), whereas the average Rpeak is 0.96, ranging 
from 0.36 (NGC 5257) to 3.21 (NGC 6240). On aver- 
age, Rpeak is two timcs higher than Rtotai- The two 
ratios are comparable for some of the compact sources 
such as IRAS 17208-0014, UGC 5101, Arp 299W. or 
NGC 5331N, but much larger in NGC 6240 (a factor of 
four larger). The lowest Rpeak is derived in the brightest 
ULIRG in our sample, IRAS 17208-0014. The general 
inference of these results is that the CO (3-2) transition 
is nearly thermalized (in the optically thick limit) in the 
inner 1-3 kpc of the U/LIRGs, whereas the CO (3-2) 
transition is subthermally populated in the extended out- 
skirts. This result is consistent with the theoretical pre- 
diction by Narayanan et al. (2007b), possibly suggesting 
that the scatter seen in the U/LIRGs in Figure 1 reflects 
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the excitation characteristics of each galaxy. FinaUy, we 
note that the total CO (3 2) emission in SMGs/quasars 
is nearly thermalized (Weiss at al. 2007b), and this could 
be the reason for the tighter correlation seen in the SMGs 
in Figure 1. 

In order to obtain a better physical understanding of 
the observed variation in a-moi as well as the observed 
scatter in the SFEs, we constructed a series of Large Ve- 
locity Gradient (LVG) models (Goldrcich & Kwan 1976) 
using the publicly available RADEX code (van der Tak 
et al. 2007). Figure 3 shows the expected flux as a func- 
tion of H2 density (sec figure caption for the adopted 
input parameters). Three key arguments can be made 
from this figure. First, the rise in flux is nearly linear up 
to ~ W cm"3 foj. CO (3-2) and up to ~ 10^ cm'^ for 
HON (1-0), suggesting that these molecular transitions 
are good tracers of the density of the molecular medium 
up to these characteristic densities. The predicted flux 
beyond these characteristic values reaches an asymp- 
totic value as the relevant lines become thermalized (i.e. 
Tex ~ T kinetic), thus predicting the gas densities from 
these lines alone becomes difficult in this regime. Sec- 
ond, the predicted flux of CO (3-2) and CO (1-0) hnes 
cross at ~ lO'' cm~^ (i.e. flux ratio of unity), and the 
line ratio is less than unity for densities below this critical 
value. Lastly, the temperature dependence of CO (3-2) 
is significant compared to the CO (1 0) line because the 
excitation temperature is higher in CO (3 2). At low 
densities (riHa ~ 10^ cm^-'), the observed CO (3 2) flux 
can vary by two orders of magnitude if the temperature of 
the ISM varies from galaxy to galaxy. Thus, at densities 
less than the critical density of the transition, coUisional 
and radiative excitation both play key roles in populating 
the energy levels. 

Under the assumptions adopted for the LVG model 
used here, the observed variation in amol can be ex- 
plained in the following way. The low Q!co(i-o) likely 
reflects the nature of CO (1-0) emission as a low density 
tracer, and thus the higher density gas in U/LIRGs do 
not further increase the observed CO (1-0) flux, leading 
to the flattening of Q!co(i-o)- The CO (3-2) emission is 
likely a tracer of higher density gas, but its flux is also 
sensitive to temperature, most notably in the low den- 
sity end. Therefore, the relatively large scatter in the 
U/LIRGs and local galaxies in Figure 1 (left) could be 
explained by intrinsic variation in the average gas kinetic 
temperature and density, reflecting the physical charac- 
teristics of each galaxy. The results from LVG modeling 
are qualitatively consistent with the theoretical predic- 
tions by Krumholz & Thompson (2007) and Narayanan 
et al. (2007b). 

Finally, we note that because the excitation of the 
CO (3-2) line coifld be sensitive to temperature, less 
temperature dependent higlic;r density gas tracers such 
as HON (1-0) or HCO+ (Krumholz & Thompson 2007; 
Gracia-Carpio et al. 2008, but see Papadopoulos (2007) 
for a counter argument for the use of HCO"'") could be 
the better tracer of star formation in quiescent galaxies, 
although the gas must be dense (tihs ~ 10^~^ cm~^) 
enough to excite the HCN emission. These lines are also 
much weaker than CO lines and therefore harder to ob- 
serve (see e.g. Figure 3). While HCN emission is often 
used a tracer of dense gas feeding star formation activity. 



HCN emission from X-ray irradiated circumnuclear tori 
may become significant and therefore the HCN emission 
may be more pronounced towards AGN than starbursts 
in some cases (Kohno et al. 2001). 

4. THE EVOLUTIONARY STATUS OF SUBMILLIMETER 

GALAXIES 

4.1. CO (3-2) Source Size 

In §3, we have shown evidence that the overall CO (3- 

2) luminosity is well correlated with the total dusty star 
formation traced in FIR luminosities. In this section, we 
extend this argument and assume that the large scale 
spatial distribution of CO (3-2) emission is also corre- 
lated with the extent of current star formation. To this 
end, we have plotted the CO (3-2) emitting diameter (in 
kpc) against the FIR luminosities in Figure 4. For the 
U/LIRGs, a trend is seen where the CO (3 2) size de- 
creases as a function of FIR luminosity. The CO (3-2) 
size, however, can vary by an order of magnitude (0.3- 
3.1 kpc) for the LIRG population alone. In contrast, the 
ULIRGs are predominantly compact (0.4-1 kpc), and the 
large FIR luminosity seen in ULIRGs is likely related to 
this high central concentration of high density gas. 

This trend immediately breaks down in SMGs where 
the derived CO (3-2) sizes (3-16 kpc) are roughly an or- 
der of magnitude larger than the ULIRGs. The average 
CO size for the SMGs (8 ± 2 kpc) is only a factor of 
two smaller than the average nuclear separation (15 ± 6 
kpc)of the pre-coalescence LIRGs (e.g. VV 114, Arp 299; 
see also Figure 4). The three resolved quasars have 
source sizes (2-4 kpc) that are systematically smaller 
than the SMGs, but the sample size is too small for 
this trend to be conclusive. We caution that, although 
the SMGs/quasars are observed with the highest angular 
resolution achievable with existing interferometers, these 
sources are only marginally resolved, and the resultant 
source sizes should be treated as a conservative upper 
limit to any internal structure present in these galaxies. 
The upper limits for the unresolved SMGs and quasars 
(see Figure 4 (right)) do not place very useful limits on 
their intrinsic sizes, but are generally consistent with the 
conclusions drawn from the resolved sources. 

These sizes show that the star forming regions in local 
ULIRGs are much more compact than those in SMGs. 
This result is an interesting contrast to the analysis by 
Bouche et al. (2007) where they find that the millime- 
ter or CO size of the SMGs themselves are much more 
compact than the optical size of the UV or optically se- 
lected population at 2; ~ 2. The widely extended (as 
opposed to nuclear) star formation in SMGs was sug- 
gested previously using high resolution radio observa- 
tions of spectroscopically identified SMGs by Chapman 
et al. (2004a) who find 70% of sources with size scale of 
7± 1 kpc. More recently, Biggs & Ivison (2007) combined 
the radio emission data of SMGs obtained using Multi- 
Element Radio Linked Interferometer Network (MER- 
LIN) and VLA in the u-v plane, and found an average 
SMC source size of 5 kpc, which is in agreement with the 
resiflts of Chapman et al. (2004a) and also the CO size 
measurements by Tacconi et al. (2006). Two-arcsecond 
resolution SMA 890 /zm measurements (lono et al. 2006b; 
Yoimger et al. 2007, 2008) of bright SMGs suggest source 
sizes of < 8 kpc, and new higher angular resolution SMA 
measurements suggest a smaller source size of 5 kpc 



7 



(Younger et al. 2008). Thus, radio, CO, and submm ob- 
servations all suggest that the size scale of the brightest 
SMGs are of order ~ 5 kpc, and this is a factor of 2 - 
10 times larger than the CO (3-2) derived source sizes 
in U/LIRGs. The average size of the SMGs in our CO 
(3-2) selected sample (8 ±2 kpc) is consistent with these 
numbers. 

4.2. Merger Scenarios for SMGs and quasars 

In Figure 5, we plot the CO (3-2) diameters as a func- 
tion of the CO (3-2) FWHM. Roughly two thirds (18/26) 
of the source sizes for SMG / quasars arc only upper lim- 
its constrained from the observations. It is evident from 
Figure 5 that oven though the size scales of SMGs are 
comparable to the nuclear separation of the widely sep- 
arated LIRGs (i.e. + sign in Figure 5), the FWHMs are 
on average 50% larger (and a factor of two larger than 
those of the ULIRGs). 

The fact that both the average CO (3 2) sizes and the 
SFEs are similar for the SMGs and the local LIRGs sug- 
gests that they may be in a similar evolutionary state, 
e.g., that many of the SMGs are in fact early or inter- 
mediate rather than late stage mergers. Both the lo- 
cal ULIRGs (which are known to be advanced mergers) 
and the quasars are more compact (significantly so for 
the ULIRGs) and have SFEs that are a factor of three 
larger than those of the SMGs. In their analysis of several 
SMGs observed with the best angular resolution, Tacconi 
et al. (2008) also conclude that SMGs are major merg- 
ers in various evolutionary stages. They discuss three 
specific SMGs (two of which, SMM J163650-F4057 and 
SMM J123707+6214, are included in our analysis) that 
each show two components separated by 8-25 kpc, which 
is more characteristic of earlier stage mergers of which 
local examples are Arp 299 and VV 114. Given that 
most SMGs are not observed at quite such high resolu- 
tion, it remains possible that many SMGs could contain 
relatively compact CO reservoirs in two individual merg- 
ing components, but not be sufficiently separated to be 
resolved with the existing data. However, this possibil- 
ity would still be consistent with our interpretation of 
the SMGs as predominantly earlier stage mergers, as the 
resulting projected nuclear separation of SMGs would 
likely be as large as the widely separated LIRGs in our 
local sample. 

In mergers of two comparable mass galaxies, the 
FWHM derived from the total CO (3-2) spectrum con- 
tains information pertaining to both galaxy rotation and 
encounter velocity. In principle, the encounter velocity 
of two comparable mass field galaxies in parabolic or- 
bit is AV ~ \/2 Vc where AV and Vc are the encounter 
velocity and circular velocity, respectively. For local sys- 
tems, AV amounts to ^ 300 km s^^ at maximum. In 
reality, the observed AV" is significantly smaller than this 
value, owing to, for example, projection and degeneracy 
in the merger phase, such as those seen in NGC 5331 
(AV = 60 km s-i) or NGC 5257/8 {AV = 11 km s'^). 
At first glance, the broad CO (3-2) Hnewidths of SMGs 
(i.e. 600 km s^^) argue against widely separated merg- 
ers. However, a simple analysis that assumes two equal 
mass field galaxies is likely inapplicable here because the 
brightest submillimeter sources are evidently embedded 
deep in the proto-cluster potential (Stevens et al. 2003; 
Greve et al. 2005). Thus, the encounter velocity of these 



galaxies could be significantly larger than those of lo- 
cal field galaxies. For comparison, a similarly large ve- 
locity dispersion is observed in the galaxy members of 
the Virgo Cluster (Rubin, Waterman & Kenney 1999). 
The CO linewidth can trace the underlying stellar as 
well as the non-baryonic halo potential (Narayanan et al. 
2008), and the suggested average stellar mass of SMGs 
(2.5 X 10^^ Mq; Borys et al. 2005) is significantly more 
massive than a typical field galaxy in the local universe 
(i.e. the Milky Way). Therefore, we suggest that the 
larger CO (3 2) FWHMs seen in SMGs are likely due to 
a combination of an underlying more massive dark mat- 
ter potential and intrinsically more massive and gas-rich 
merger progenitors. 

Analysis of numerical simulations (Greve & Sommer- 
Larsen 2008; Narayanan et al. 2006, 2007a, 2008) have 
found that the CO profile can display a variety of shapes 
owing to strong central outfiows which could appear as 
secondary emission peaks in the CO spectra. These sim- 
ulations further predict that the spectral characteristics 
have strong dependences on the sightline (i.e. the incli- 
nation of the source) , as well as the mass of the under- 
lying potential especially toward the latter stage of the 
evolution when the gas has settled in a disk (i.e. quasar 
phase). In addition, the simulated CO profile evolves as a 
function of merger age, where kinematic asymmetries are 
seen primarily in the pre-coalcscence phase when the two 
galaxies have not yet merged. In contrast, the spectrum 
becomes closer to Gaussian as the molecular gas becomes 
virialized and a quiescent molecular disk is formed. 

While some of their adopted assumptions may not be 
best suited for analyzing the observational data pre- 
sented in this study, these simulation results provide 
important physical insight. The narrower CO (3-2) 
Hnewidths, and the dominance of Gaussian shapes in the 
quasar line profiles are consistent with the characteris- 
tics of a post-coalescence galaxy seen in simulations. In 
our limited sample of only three quasars, the size scale of 
quasars with robust size estimates are intermediate be- 
tween the nuclear separation of pre-coalescence LIRGs 
and the very compact ULIRGs (see Figure 5). Given the 
factor of 10 larger average LpiR in quasars compared to 
the local ULIRGs, we might naively expect the starburst 
region to be larger in the quasars, which could explain the 
size difference to within a factor of two. Alternatively, 
the quasar activity can begin in one of the progenitor 
galaxies earlier in the merger evolution than that pre- 
dicted in simulations, and the abmidance of CO (3 2) 
c;niitting molecular gas suggests intense concurrent star- 
burst activity. It is also possible that CO (3-2) emission 
entrained in powerful AGN jets can lead to an apparent 
increase in the CO (3-2) size of quasars (Narayanan et 
al. 2006). Further, the observed narrower CO hnewidths 
of quasars may simply be due to an intrinsically less mas- 
sive quasar disk (Greve et al. 2005) or due to geometri- 
cal effects (CarilU & Wang 2006; Wu 2007) with quasars 
preferentially discovered in more face-on galaxies. It is, 
however, impossible to obtain a robust conclusion from 
our current compilation of data which is limited in sam- 
ple size, sensitivity, and angular resolution. 

Results from numerical simulations that include nega- 
tive feedback from starbursts as well as the central AGN 
(Hopkins et al. 2005) suggest that SMGs and quasars are 
not in the same evolutionary stage, but linked through 
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dynamical evolution where the quasar activity becomes 
visible during the final stages of the coalescence. We 
have provided evidence in this paper that the CO (3-2) 
properties of many of the SMGs are consistent with the 
widely separated LIRGs in the local universe. A possible 
scenario suggested from these results is that SMGs are 
widely separated gas-rich mergers taking place near the 
centers of massive halos. Global disk-wide starburst ac- 
tivity in gas-rich high redshift mergers is also predicted 
theoretically by Mihos (2001). On the other hand, some 
of the CO (3-2) derived properties in quasars, such as 
size and SFE, are more consistent with ULIRGs than 
with LIRGs or SMGs. Overall, the data summarized 
here are broadly consistent with this evolution scenario. 

Finally, we note that there is a selection effect in play. 
There are interacting/merging galaxies in the local uni- 
verse that are not U/LIRGs; thus, the triggering of the 
U/LIRG activity also depends on the amount of pre- 
existing gas available in the merging pair. We also note 
that there are ULIRGs that are in the early stage of the 
interaction (Dinh-V-Trung et al. 2001), as well as LIRGs 
that are apparently isolated (i.e. NGC 1068). There is 
also a case like the z = 4.7 quasar - SMG pair BR1202- 
0725, in which two submm bright sources (one of which 
is a quasar, while the companion appears to harbor an 
obscured AGN, lono et al. 2006c) are separated by a pro- 
jected distance of 26 kpc. The true evolution, therefore, 
is much more complicated than we envision in these sim- 
ple scenarios, and sensitive high resolution ALMA sur- 
veys of SMGs, quasars, as well as quiescent high redshift 
galaxies in a variety of evolutionary stages should reveal 
their true nature. Lastly, the gas consumption time for 
SMGs using the currently available fuel is estimated to 
be ~ 40 Myr (Grcve et al. 2005), and this is significantly 
smaller than the typical major merger timescale of 1 Gyr. 
If the SMGs and quasars indeed represent some stage 
of the evolution of a merger with sustained starbursts, 
then the dense gas must be replenished continuously in 
order to maintain a continuous burst of star formation 
throughout the merger evolution. 

5. SUMMARY 

In this paper, we have presented a detailed compari- 
son of the CO (3-2) emitting molecular gas between a 
local sample of luminous infrared galaxies (U/LIRGs) 
and a high redshift sample that comprises submm se- 
lected galaxies (SMGs), quasars, and two LBGs. The 
data for the local sample come from our recent Submil- 
limeter Array survey of CO (3-2) emission in U/LIRGs 
while the CO (3 2) data for the high redshift population 
are obtained from the literature. 

(1) We find that Lco(3-2) ^'^^ LpiR are correlated over 
five orders of magnitude, which suggests that the molec- 
ular gas traced in CO (3-2) emission is a robust tracer 
of star formation activity. The slope derived from the 
log Lco(3-2) ~ ^*^S LpiR relation is 0.93 ± 0.03, and this 
near unity slope suggests that the efficiency of convert- 
ing CO (3 2) emitting molecular gas to massive stars is, 
to within a factor of two, nearly uniform across different 
types of galaxies residing in different epochs. Within the 
local sample, the global star formation efficiency is lower 
in early to intermediate stage mergers than in the centers 
of more advanced stage mergers. SMGs show lower star 



formation efficiencies that are comparable to those of the 
widely separated LIRGs. 

(2) The slope derived for the Lco(3-2) ~ ^fir relation 
is significantly steeper than the slope derived from the 
same relation using the CO (1-0) emission. We show 
that this slope approaches unity for higher critical den- 
sity tracers. Further, we show that the CO (3-2) emis- 
sion in the central regions of U/LIRGs is, on average, 
nearly thermalized, whereas the outskirts are subther- 
mally populated. 

(3) From non-LTE LVG models as well as published 
merger evolution models that include radiative transfer, 
we argue that the CO (3-2) line is a fairly good tracer of 
star formation in SMGs, where the star forming gas den- 
sity is high enough that the dependence on temperature 
is relatively low. In contrast, Lco(3-2) ^^"^ show signifi- 
cant scatter when the average gas density of the medium 
is lower than the critical density, which appears to occur 
more often in widely separated mergers in our sample. 
The scatter in the Lj^q^.^ 2) ~ Lfir relation, therefore, 
reflects the characteristics of the galaxy (temperature, 
density, and AGN heating of dust in some cases), but 
over a large range of luminosities, the CO (3-2) line ap- 
pears to be a good probe of star formation. 

(4) The CO (3-2) derived source sizes in U/LIRGs 
show an apparent trend that the brighter ULIRGs are 
systematically more compact than the less FIR bright 
LIRGs. The CO (3-2) sizes of the SMGs arc on av- 
erage an order of magnitude larger than the ULIRGs, 
and are comparable to the separation of the widely sep- 
arated LIRGs in our sample. The three quasars with 
robust size constraints have source sizes that are some- 
what smaller than the SMGs and a factor of two larger 
than the ULIRGs. 

(5) The similarity in the CO (3-2) size and star forma- 
tion efficiency between SMGs and LIRGs suggests that 
many of the SMGs studied here could be intermediate 
stage mergers. The SMG linewidths are, on average, 
much broader than many of the U/LIRGs, and we ar- 
gue that the large encounter velocity likely arises from 
the massive halo potential in a proto-cluster environ- 
ment. In contrast, quasars have smaller CO (3-2) sizes 
and linewidths, and a dominance of Gaussian shapes in 
the line profiles. These characteristics are consistent with 
those seen in post-coalescence galaxies in encounter sim- 
ulations. 
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LEG 


2.727 


31.8 


10.99 


22.85 




37 


J213512-010143 


LEG 


3.074 


8 


11.48 


< 24.01 




38 


Average^ 
















U/LIRGs 




129 Mpc 




11.74 


23.37 


49 




ULIRGs 




178 Mpc 




12.26 


23.77 


61 




LIRGs 




116 Mpc 




11.60 


23.27 


45 




SMGs 




2 = 2.493 


5.3 


12.76 


24.19 


7.8 




Quasars 




z = 2.782 


12.8 


13.09 


24.42 


9.0 




LEGs 




z = 2.898 


19.9 


11.23 


22.85 







Note. — See Appendix of Paper I for the discription of individual sources. 

^Values for U/LIRGs are given in Dz, (Mpc) and the high redshift population are given in redshifts. For those sources that are resolved into two 
or more galaxies, the projected nuclear separation derived from 2MASS K-band images are shown in () in units of kpc. 

^Values are corrected for lensing if the lensing factor is known. See text for an explanation of how the values are derived. The FIR luminosity of 
cB58, and J213512 were obtained from Baker ct al. (2004) and Coppin et al. (2007), respectively. 

^ The 1.4 GHz data were obtained from the NVSS (Condon ct al. 1998) for the U/LIRGs, Small ct al. (2002) Chapman et al. (2005) for the SMGs, 
Carllli et al. (2004) White et al. (1997) Yun et al. (2000), Barvainis, Lonsdale & Antonueci (1996), for the quasars, and Becker et al. (1995) for the 
LBGs. Values are corrected for lensing if the lensing factor is known. 

*^SMA 880/im flux for the U/LIRGs (paper I). SCUBA 850/im flux for the high redshift galaxies. Values are corrected for lensing magnification 
when the magnification factor is known. 

^1. Paper I, 2. Frayer et al. (1998), 3. Gcnzel et al. (2003), 4. Small, Ivison, & Blain (1997), 5. Small et al. (2002), 6. Tacconi et al. (2006), 
7. Neri et al. (2003), 8. Chapman et al. (2005), 9. Frayer et al. (1999), 10. Small et al. (1998), 11. Swinbank et al. (2004), 12. Motohara et al. 
(2005), 13. Small, Smith & Ivison (2005), 14. Sheth et al. (2004), 15. Kneib et al. (2004), 16. Scott et al. (2002), 17. Greve et al. (2005), 18. 
Ivison et al. (2002), 19. Greve et al. (2004), 20. Barger, Cowie & Sanders (1999), 21. lono et al. (2006a), 22. Borys et al. (2006), 23. Walter et al. 
(2004), 24. Robson et al. (2004), 25. Mainline et al. (2004), 26. Barvainis & Ivison (2002), 27. Guilloteau et al. (1999), 28. Barvainis et al. (1998), 
29. Solomon & vanden Bout (2005), 30. Priddey et al. (2003), 31. Beelen et al. (2004), 32. Weiss et al. (2003), 33. Downes, Solomon & Radford 
(1995), 34. Rowan-Robinson et al. (1993), 35. Knudsen, van der Werf, & Jaffe (2003), 36. Coppin et al. (2008), 37. Baker et al. (2004), 38. Coppin 
et al. (2007), 39. Knudsen et al. (2008) 

^ Scaled the 1.2 mm flux in Coppin et al. (2008) to the 850/im flux by assuming /3 = 1.5. 

^ Excluding limits. 
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TABLE 2 

CO (3-2) Derived Properties 



Source 



'^CO(3-2)t 



A FWHM 



Source Size'*! 



(7^ 



FIR 



CO(3-2) 



U/LIRGs 

IRAS 17208-0014 
Mrk 231 
Mrk 273 

IRAS 10565+2448 
UGC 5101 
Arp 299 (total)'^ 

Arp 299E 

Arp 299W 
Arp 55 (total)"^ 

Arp 55N 

Arp 55S 
Arp 193 
NGC 6240 
VV 114 

NGC 5331 (totaI)<= 

NGC 5331S 

NGC 5331N 
NGC 2623 
NGC 5257/8 (total)'^ 

NGC 5257 

NGC 5258 
NGC 1614 

High Redshift Galaxies 

SMMJ02399-0136 
SMMJ0443007+0210 

SMMJ123549+6215 

SMMJ123707+6214 (total) 
SMMJ123707+6214SW 
SMMJ123707+6214NE 

SMMJ14011+0252 



h6612 
-4057 
-4105 
h4053 
[-0015 



-661144 



SMMJ16359-I 
SMMJ163650H 
SMMJ163658H 
SMMJ1637H 
SMMJ22174-I 
MIPS-J1428 
SMMJ16354H 
SDSSJ1148+5251 
RXJ0911. 4-1-0551 
LBQS1230-I-1627B 
MG0414-I-0534 
LBQS0018 
VCVJ1409+5628 
Clovcrleaf 
IRASF10214-I-4724 
SMMJ04135-I-1027 
HS1002-I-4400 
RXJ124913-055906 
SMMJ131444+423814 
VV96 J140955-I-562827 
VV96 J154359-I-535903 
HS1611+4719 
MS1512-CB58 
J213512-010143 

Average' 

U/LIRGs 

ULIRGs 

LIRGs 
SMGs 
Quasars 
LBGs 



4.5 X 10^ 

2.6 X 10^ 

3.2 X 10^ 
1.9 X lO'' 

1.7 X 10^ 

1.3 X 10^ 

8.4 X 10* 
3.1 X 10* 
1.1 X 10^ 
7.4 X 10* 
4.1 X 10* 
2.4 X 10^ 
7.4 X 10^ 

3.1 X 10^ 

2.6 X 10^ 

2.2 X 10^ 

4.3 X 10* 
1.0 X 10^ 
3.0 X 10^ 

6.7 X 10* 
2.2 X 10^ 

8.6 X 10* 

4.7 X IQio 
1.0 X IQio 
4.2 X IQio 
2.9 X 10^° 

1.9 X 101° 
1.0 X 10^0 
1.9 X 10^° 
3.6 X lO^s 
6.9 X lOio 

5.6 X 10^° 

3.0 X 10^° 

3.7 X lOio 

1.4 X 10" 

2.8 X 101° 

2.6 X 10^° 

5.1 X 10^ 
3.0 X 101° 

9.2 X 10i° 
5.4 X 10i° 

7.9 X 101° 

4.0 X 10i° 

6.7 X lO'' 
1.7 X lOii 
4.2 X 10i° 
3.6 X 10i° 
1.2 X 10i° 
8.2 X 10i° 

3.1 X 10i° 
5.1 X 10i° 
4.4 X 10* 
2.9 X 10^ 



(2.6 ±0.5) X 10^ 
(3.4 ± 0.6) X 10^ 
(2.4 ± 0.6) X 10^ 
(4.4 ±1.1) X lOio 
(5.0 ± 1.0) X 10i° 
(1.6 ± 1.2) X 10^ 



570 
200 
500 
200 
650 

400 
250 

200 
200 
380 
500 
300 

450 
370 
400 

250 

370= 
280 



1360 ± 50 
350 ± 40 
600 ± 50 

430 ± 60 
430 ± 60 
190 ± 11 
500 ± 100 
710 ± 50 
800 ± 50 
830 ± 130 
780 ± 100 
386 ± 104 
284 ± 50 

280 
350 ± 60 

580 
163 ± 29 
311 ±28 
416 ± 6 
220 ± 30 
340 ± 120 
640 ± 160 
1090 ± 340 
550 ± 220 
310 ±30 
520 ± 140 
230 ± 40 
174 ± 43 
190 ± 24 



360 ± 30 
420 ±110 
350 ± 30 
590 ± 90 
430 ± 60 
180 ± 10 



< l.Ox < 0.8 
0.5 X 0.3 
0.4 X 0.3 

1.0 X 0.7 

1.1 X 0.8 

0.3 X 0.2 
0.4 X 0.3 

0.5 X 0.5 

< 0.8x < 0.5 
1.6 X 0.4 
0.9 X 0.6 

3.0 X 2.2 

1.5 X 0.9 

< 1.9x < 0.8 
0.5 X 0.5 

< 3.1X < 1.6^ 

3.1 X 2.6<i 
1.0 X 0.8 
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< 3.1X < 0.7 
<4.2 (2.5)f 

7.3 

< 4.0 
11 

12x < 6 
6.5x < 3.3 
3.2 

< 60x < 41 

< 69x < 50 

< 36 

< 57x < 41 
3.6x < 1.4 

< 22x < 12 

< 48x < 24 

< 24 

< 20x < 13 

< 20x < 12(2.4) 
3.8x < 2.5 

< 96x < 48 

< 38x < 28 

< 73x < 54 

< 42x < 33 

< 19x < 14 

< 46x < 29 

< 41x < 33 

< 65x < 38 (2.0)'' 

< 48x < 42 



(1.1 X 0.8) ±0.3 
(0.5 X 0.3) ±0.1 
(1.3 X 0.9) ±0.3 
8.4 ± 1.9 
3.3 ±0.4 



564 (2.75) 
793 (2.90) 
378 (2.58) 
439 (2.64) 
448 (2.65) 
409 (2.61) 



349 (2.54) 



159 (2.20) 

47 (1.68) 
102 (2.01) 
117 (2.07) 



291 (2.46) 
93 (1.97) 



314 (2.50) 

203 (2.31) 
171 (2.23) 
149 (2.17) 
206 (2.31) 



164 (2.21) 
196 (2.29) 
113 (2.05) 
110 (2.04) 
262 (2.42) 
166 (2.22) 
173 (2.24) 
536 (2.73) 
1000 (3.00) 
222 (2.35) 
816 (2.91) 
168 (2.23) 
227 (2.36) 
1084 (3.03) 
188 (2.27) 
804 (2.91) 
351 (2.55) 
262 (2.42) 
194 (2.29) 
250 (2.40) 
329 (2.52) 
323 (2.51) 
235 (2.37) 
224 (2.35) 
105 (2.02) 



322 ± 56 
578 ± 120 
252 ± 46 
194 ± 20 
430 ± 83 
170 ± 60 



Note. — Units are [K km s ^ pc^] for log Lco(3 — 2)' s ^] for A FWHM, [kpc] for source size, and [L0 (K km s"^ pc-^)"^)] for 

Lfir/Lco(3-2) linear and logarithmic scales. 

^ The deconvolved source size. Gaussian fits for IRAS 17208-0014, Arp 55S, NGC 5331N, NGC 5257 give only the upper limits shown. Sizes 
derived from CO transitions other than CO (3—2) are used for some of the high redshift sources. For high redshift sources with known luminosity 
magnification factors, we have scaled the CO size by the square root of the luminosity magnification assuming that the surface brightness is conserved 
during gravitational lensing. 
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TABLE 3 

CO (3-2) - CO (1-0) Ratios in U/LIRGs 





p a 
^total 


R t> 
^peak 


CO (1-0) 


IRAS 17208-0014 


0.40 


0.41 


1 


Mrk 231 


0.35 


0.82 


1 


Mrk 273 


0.63 


... d 


1 


IRAS 10565+2448 


0.33 


0.72 


1 


UGC 5101 


0.46 


0.48 


2 


Arp 299E 


0.74 


1.11 


3 


Arp 299W 


1.25 


1.30 


3 


Arp 55N 


0.20 


... 6 


4 


Arp 55S 


... e 


4 


Arp 193 


0.61 


< 1.26 


1 


NGC 6240 


0.83 


3.21 


5 


VV 114 


0.25 


0.49 


6 


NGC 5331S 


0.30 


0.55 


7 


NGC 5331N 


0.50 


0.50 


7 


NGC 2623 


0.44 


0.74 


5 


NGC 5257 


0.21 


0.36 


7 


NGC 5258 


0.38 


1.06 


7 


NGC 1614 


0.35 


f 


8 


average 


0.48 


0.93 





^ The CO (3-2) to CO (1-0) flux ratio using the total luminosities. 

^ The CO (3-2) to CO (1-0) flux ratios using the peak fluxes. The peak fluxes were derived after convolving the CO (3-2) maps to the lower 
resolution CO (1-0) maps, except for Arp 193 where the CO (1—0) map has higher angular resolution than the CO (3—2) map. The peak CO (1—0) 
fluxes were obtained from the images published in the references. These ratios trace the inner 1—3 kpc of the galaxies. 

*^ 1. Downes & Solomon (1998), 2. Gcnzel ct al. (1998), 3. Aalto et al. (1997) 4. Sanders et al. (1988b), 5. Bryant &: Scoville (1999), 6. Yun et 

al. (1994) 7. lono, Yun & Ho (2005), 8. Scoville ct al. (1989) 

*^ There appear to be inconsistencies in the published values of Mrk 273 in Fig. 15 of Downes & Solomon (1998). The peak flux level (~ 
144 Jy km s""'" beam""'") exceeds the total flux (78 Jy km s""*") given in the same paper. 

^ Although the CO (1—0) contours arc shown overlaid on an optical image, the contour levels were not explicitly given in Sanders et al. (1988b). 

^ NGC 1614 is unresolved in this paper. 



